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Abstract 

We study a model for itinerant, strongly interacting fermions where a judicious tuning of the 
interactions leads to a supersymmetric Hamiltonian. On the triangular lattice this model is known 
to exhibit a property called superfrustration, which is characterized by an extensive ground state 
entropy. Using a combination of numerical and analytical methods we study various ladder ge- 
ometries obtained by imposing doubly periodic boundary conditions on the triangular lattice. We 
compare our results to various bounds on the ground state degeneracy obtained in the literature. 
For all systems we find that the number of ground states grows exponentially with system size. 
For two of the models that we study we obtain the exact number of ground states by solving the 
cohomology problem. For one of these, we find that via a sequence of mappings the entire spectrum 
can be understood. It exhibits a gapped phase at 1/4 filling and a gapless phase at 1/6 filling and 
phase separation at intermediate fillings. The gapless phase separates into an exponential number 
of sectors, where the continuum limit of each sector is described by a superconformal field theory. 
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I. INTRODUCTION 



The supersymmetric model for lattice fermions was first introduced in [T3] and has sub- 
sequently seen an intensive follow up in both the physics and mathematics literature (for 
a fairly recent overview see |18| and references therein). The key property of this model, 
namely supersymmetry, has the remarkable consequence that exact results can be obtained 
in the non-perturbative regime of a model for strongly interacting itinerant fermions. 

These exact results have revealed a wide range of interesting phenomena in this system. 
On one dimensional lattices one can typically find quantum criticality, which in combination 
with supersymmetry, gives rise to an effective continuum limit description in terms of a 
superconformal field theory [5l HSf [151 HH]- On general lattices the system furthermore 
exhibits a strong form of quantum charge frustration called superfrustration jTH |22| l33] . 
This property is characterized by an extensive ground state entropy. Other phenomena, such 
as edge modes, possible topological order and Rokhsar-Kivelson like quantum liquid phases 
were also observed for this model fIU\ fH\ 123] . Finally, a surprising connection between the 
supersymmetric model and the XYZ spin chain along a certain line of couplings has made it 
possible using supersymmetry to reinterpret some known properties and derive new results 
for the spin chain [Tn[T2l[T7]. 

On the Mathematics side the study of supersymmetric lattice models has led to interesting 
results on the cohomology of independence complexes of lattices and graphs, 2D grids in 
particular [21 |U [HI [21H2S]- The two sides are connected by the observation that 
quantum ground states of the supersymmetric lattice model are in 1-to-l correspondence 
with the elements of the cohomology of an associated independence complex. 

In this paper we focus on the supersymmetric model on the triangular lattice. The 
ground state structure of this model is not fully understood. Various techniques that were 
successfully used to explore the ground state structure on other lattices have failed so far 
for this lattice. Nevertheless, a number of non-trivial results have been established for this 
model. The numerical computation of a lower bound on the number of ground states has 
shown that the system has an extensive ground state entropy |33]. Furthermore, analytic 
results on cohomology elements have revealed that degenerate quantum ground states occur 
in the entire range between 1/7 and 1/5 filling [28j. Finally, an upper bound on the number 
of quantum ground states was also obtained analytically P|. 
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In this work we focus on ladder realizations of the triangular lattice by imposing dou- 



bly periodic boundary conditions (see figure 1. 1 ) . Using a combination of numerical and 
analytical techniques we study the ground state structure of these systems. We compare 
our results to the bounds that were previously established. For all systems we find that 
the number of ground states grows exponentially with system size. We numerically confirm 
that the range of filling for which zero energy states exist agrees with the range obtained 
by Jonsson. Furthermore, we prove analytically that there are no zero energy states outside 
the range between 1/8 and 1/4 filling. For two of the models that we study we obtain the 
exact number of ground states by solving the cohomology problem. 

For the simplest ladder, the 2 leg ladder, we uncover a symmetry that allows us to 
understand the origin of the exponential ground state degeneracy. This local Z2 symmetry 
distinguishes an exponential number of sectors. Subsequently, we present a mapping from 
the ladder to the chain such that the entire spectrum can be understood. Remarkably, in 
certain sectors the spectrum turns out to be gapped, in others it is gapless and in yet other 
sectors we find phase separation. The continuum limit of each of the gapless sectors is shown 
to be described by a superconformal field theory with central charge c = 1. 

The paper is organized as follows. In the remainder of this section, we define the super- 
symmetric model, discuss the analytic and numerical techniques that we employ and, finally, 
summarize the results on the ground state degeneracy of the supersymmetric model on the 
triangular lattice that were established previously. In section |n| we present the numerical 
results we obtained for systems up to 54 sites and discuss how they relate to the known 



bounds and what they imply for the full 2D triangular lattice. In section |III[ we focus on 
the 3 and 4 leg ladders for which we obtain a number of analytical results on the ground 



state structure. Finally, in section IV we present the full solution of the 2 leg ladder. 



A. The model 



Supersymmetry is a symmetry between fermionic and bosonic degrees of freedom (see |3] 
for a general reference). It plays an important role in theoretical high energy physics, where 
various theories that go beyond the standard model require supersymmetry for a consistent 
formulation. In these theories all the known elementary particles are accompanied by yet to 
be discovered superpartners. 



3 



Figure I.l: In this paper we consider the supersymmetric model on the triangular lattice with 
doubly periodic boundary conditions. For the M x N triangular lattice the boundary conditions 
are imposed by requiring periodicity along two directions: u = (0,-/Vf) and v = {N,0). The figure 
shows the two unit vectors along which the periodicities are imposed. 

In the lattice model discussed here the physical particles are spinless lattice fermions and 
the supersymmetry relates fermionic and bosonic many particle states with an odd and even 
number of the lattice fermions, respectively. In the M = 2 supersymmetric theories that 
we consider, a central role is played by the operators Q and Q\ called supercharges, which 
have the following properties [33] 

• Q adds one fermion to the system and takes out one fermion from the system. 

• The supercharges are fermionic operators and thus nilpotent: = (Q^)^ = 0. 

• The hamiltonian is the anti-commutator of the supercharges, H = {Q,^^}, and as a 
consequence it commutes with the supercharges and conserves the number of fermions 
in the system. 

Imposing this structure has some immediate consequences: supersymmetric theories are 
characterized by a positive definite energy spectrum and a twofold degeneracy of each non- 
zero energy level. The two states with the same energy are called superpartners and are 
related by the supercharge. 

From its definition it follows directly that H is positive definite: 

= iQiv^)r+iQ^i^)r>o. 

Furthermore, both Q and commute with the Hamiltonian, which gives rise to the 
twofold degeneracy in the energy spectrum. In other words, all eigenstates with an energy 
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Eg > form doublet representations of the supersymmetry algebra. A doublet consists of 
two states {\s),Q\s)), such that Q^s) = 0. Finally, all states with zero energy must be 
singlets: Q\g) = Q^g) = and conversely, all singlets must be zero energy states [34]. In 
addition to supersymmetry our models also have a fermion-number symmetry generated by 
the operator F with 

[F,Qt] = _gt and [F,Q]=Q. (I.l) 

Consequently, F commutes with the Hamiltonian. Furthermore, this tells us that superpart- 
ners differ in their fermion number by one (let F\s) = fs\s), then F{Q\s)) = Q{F + l)\s) = 

(/, + i)(g|.))). 

An important issue is whether or not supersymmetric ground states at zero energy occur, 
that is, whether there are singlet representations of the algebra. For this one considers the 
Witten index 

W = tT[{-lfe-^''] , (1.2) 

where the trace is over the entire Hilbert space. Remember that all excited states come in 
doublets with the same energy and differing in their fermion-number by one. This means 
that in the trace all contributions of excited states will cancel pairwise, and that the only 
states contributing are the zero energy ground states. We can thus evaluate W in the limit 
of /3 — )■ 0, where all states contribute (—1)^. It also follows that \W\ is a lower bound to 
the number of zero energy ground states. 

We now make things concrete and define a supersymmetric model for spin-less fermions 
on a lattice, following [T^. The operator that creates a fermion on site i is written as c| with 
{c\,Cj} = 6ij. To obtain a non-trivial Hamiltonian, we dress the fermion with a projection 
operator: P<j> = Ylj next to «(-'- ~ which requires all sites adjacent to site i to be empty. 
With Q = '^iP<i> = X] CiP<^iy, the Hamiltonian of these hard-core fermions reads 

H = {QlQ} = Y^ E P<i>4c,P<J> + Y.P<^>. (1.3) 

i j next to i i 

The first term is just a nearest neighbor hopping term for hard-core fermions, the second 
term contains a next-nearest neighbor repulsion, a chemical potential and a constant. The 
details of the latter terms will depend on the lattice we choose. In particular, for the 
triangular lattice we find 

5^P<,>=Ar-6F + ^\/(,), (1.4) 

i i 
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where V(^i) + 1 is the number of particles adjacent to site i. 

B. Methods 

1. Cohomology techniques 

For the supersymmetric models, cohomology has proven to be a very powerful tool to 
extract information about the zero energy ground state (s) of the models (see for example 
[T^ [T5l [T8l [20l EH EH] and references therein). The key ingredient is the fact that ground 
states are singlets; they are annihilated both by Q and Q^. This means that a ground state 
l^f) is in the kernel of Q: Q\g) = and not in the image of Q, because if we could write 
Ifl') = Ql/); then {\f),\g)), would be a doublet. So the ground states span a subspace Hq of 
the Hilbert space H of states, such that Hq = keiQ/lmQ. This is precisely the definition 
of the cohomology oi Q. So the ground states of a supersymmetric theory are in one-to-one 
correspondence with the cohomology of Q. It follows that the solution of the cohomology 
problem gives the number of zero energy states for each particle number sector. Equivalently, 
we find that zero energy states are in one-to-one correspondence with the homology elements 

of gt. 

We compute the cohomology using the 'tic-tac-toe' lemma of [H]- This says that under 
certain conditions, the cohomology Hq for Q = Qi + Q2 is the same as the cohomology of 
Qi acting on the cohomology of Q2. In an equation, T-Lq = T-LQ^^iT-LQ^) = 7^12, where Qi and 
Q2 act on different sublattices Si and 6*2. We find by first fixing the configuration on 
all sites of the sublattice 5*1, and computing the cohomology Hqj. Then one computes the 
cohomology of Qi, acting not on the full space of states, but only on the classes in "Hqj- ^ 
sufficient condition for the lemma to hold is that all non-trivial elements of ?^i2 have the 
same /2 (the fermion-number on 6*2). 

Although ground states are in one-to-one correspondence with cohomology elements, the 
two are not equal unless the cohomology element happens to be a harmonic representative 
of the cohomology. Harmonic representatives are elements of both the cohomology of Q and 
the homology of Q^. So they are annihilated by both supercharges, which is precisely the 
property of a zero energy state. It follows that, although the solution of the cohomology 
problem gives the number of ground states, it typically does not give the ground states 
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themselves. For a more leisurely introduction to cohomology and an exposition of the 
relation between the supersymetric model and independence complexes we refer the reader 
to [BEl]. 

We now briefly state the cohomology results for the one dimensional chain [TH], because 
we will use these results many times throughout the paper. For the periodic chain of length 
L = 3?2 + a the cohomology is trivial for all particle numbers, /, except for f = n, where we 
have 

{1 for a = ±1, 
(1.5) 
2 for a = 0. 

Similarly, for open boundary conditions we have 

{1 for a = and a = —1, 
(1.6) 
for a = 1, 

for a chain of length L = 3n + a and n particles and dim('HQ) = at all other particle 
numbers. In particular, we note that the cohomology of an isolated site that can be both 
empty and occupied (L = 1 with open boundary conditions) is trivial. This is equivalent 
to the statement that the single site chain has no zero energy states, indeed the empty and 
occupied state form a doublet of energy E = 1. 

2. Numerical methods 

To complement the analytical techniques we also investigated a range of numerical meth- 
ods. The goal of these techniques is to find accurately the number of degenerate zero 
energy eigenvalues of a given Hamiltonian. In principle there are many ways to compute 
this number. Here we describe some algebraic approaches we tried and also the iterative 
diagonalization method which was found to be the most effective. 

1. Characteristic polynomial: Since the roots of the characteristic polynomial of a given 
matrix are the eigenvalues of the matrix, it is instructive to obtain the characteristic polyno- 
mial equation of the above Hamiltonian and solve it. However, the characteristic polynomial 
of a matrix bigger than 4 x 4 is a univariate of degree more than 4 which is known to be 
exactly non-solvable in terms of the radicals of its coefficients as the Abel-Ruffini theorem 
states. This means that we have to rely on numerical methods to solve univariate polyno- 
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mials for the Hamiltonians of sizes bigger than 4x4. Such numerical methods, e.g., the 
Newton's method, are impractical because even small round-off errors in the coefficients of 
the characteristic polynomial can end up being a large error in the eigenvalues and hence in 
the eigenvectors. 

However, since in this specific problem we know that the ground state eigenvalue is 
exactly zero and we thus only need to compute its degeneracy, we may get away from the 
round-off error problem by only obtaining the characteristic polynomial and studying its 
structure. In particular, since we know that the lowest lying eigenvalues in our case are 
all degenerate eigenvalues, we can do the following: obtain the characteristic polynomial of 
the given Hamiltonian of size, say, A^. Then, factorize the univariate polynomial (in the 
variable, say, x) such that the final form of the polynomial is x"^{akX^ + ... -|- aox^) ■ ■ ■ = 0, 
where m + k = N. Here, m is now the number of the degenerate lowest lying eigenvalues. 
The problem here is that obtaining the characteristic polynomial itself is a quite extensive 
task as the matrix size increases. In particular, we could not easily go beyond the 4x6 
triangular lattice (where the largest matrix dimension is 1188) on a single computer. 

Looking at the structure of the characteristic polynomial, we can also deduce that we 
may just need the smallest i in the characteristic polynomial written in the generic form 
as ^il^OjX* for which 7^ since this i is m in the above notation, i.e., the number of 
degenerate lowest lying eigenvalues. However, again, since is the sum over all the i x i 
minors of the matrix, the computation immediately blows up restricting us again to the 
matrices of sizes as mentioned above. 

2. Rank of the matrix: the rank r of a matrix of dimension is less than if there are 
N—r zero modes. Since our Hamiltonian is known to have zero eigenvalues as the degenerate 
lowest lying eigenvalues, — r is the number of lowest lying eigenvalues in this case. This 
means that we just need to compute the rank of our Hamiltonian. However, the computation 
of the rank of a matrix is an extremely demanding calculation. We tried the singular 
value decomposition (the in-built routine of Matlab) and also computing the rank directly 
using various methods including the recently developed so-called numerical rank-revealing 
method (with the RankRev package) [29^ ^Oj . The latter method did perform slightly better 
compared to the singular value decomposition routine from Matlab. The method indeed 
correctly computes the number of degenerate eigenvalues, but scales exponentially with 
increasing matrix size and hence is not effective for even moderately sized lattices. 
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3. Row-Echelon Form: Another way of computing the number of degenerate eigenvalues 
of a given matrix is to transform the matrix into Row-Echelon form and then to count the 
number of zero rows it has. Again, obtaining the Row-Echelon form of a matrix becomes a 
very difficult task since the number of steps in the algorithm increases exponentially with 
the dimension of the matrix. In particular, for the 4x4 lattice (matrix dimension of 96) 
this method completes successfully instantaneously (using Matlab's rref command), whereas 
already for the 4x6 lattice (matrix dimension 1188) this takes hours. 

4. Iterative diagonalization: Here, we used the Krylov-Schur algorithm implementation 
from the SLEPc package [I] to compute the lowest lying eigenvalues. The DoQO (Di- 
agonalisation of Quantum Observables) toolkit was used to build the matrices and call the 
SLEPc routines. This method was chosen for its ability to resolve large amounts of de- 
generate states which is not possible with the standard Lanczos algorithm. These methods 
scale linearly with the dimension of the matrix, but exponentially with the dimension of the 
degenerate subspace one is trying to resolve. In addition, unlike the techniques described 
above, it is not always guaranteed that all the degenerate zero energy states will be found. 
In particular if the gap to the next excited state is too small or the dimension of the subspace 
within which the algorithm works is too small then the algorithm will not be able pick up 
all the zero energy degenerate states. The subspace can be increased but at the expense of 
exponentially more computational resources. 

In our calculations we have taken steps to overcome these problems. Firstly we exploit the 
symmetries of the Hamiltonian including particle conservation and space group symmetries 
to drastically reduce the dimension of the matrices involved and also the dimension of the 
degenerate zero energy subspace for each matrix. As well as this we check to ensure that the 
Witten index we compute matches that calculated via transfer matrix methods. This offers 
a very good check to ensure that we accurately account for all the zero energy ground states. 
Using these techniques we have been able to successfully calculate the numbers of degenerate 
ground states for triangular lattices with up to 54 sites (see appendix |B| corresponding to 
the 6x9 triangular lattice (the full Hilbert space dimension here is on the order of 100 million. 
With filling and momentum conservation the largest matrix is on the order of 50 thousand). 

In the methods 1 — 3, we do not need to know the Witten index at all. But there the 
computation becomes heavy. The methods 1 and 3 can be used for the symbolic matrices 
too. So if there are some parameters involved in the Hamiltonian, for example, for the 
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staggered matrices, we can still find the row-echelon form or the characteristic polynomial 
of the corresponding matrices. This is not the case for the purely numerical methods 2 and 
4. 

All in all, in our experiments with different methods, we conclude that the iterative 
diagonalization method scales better than any of the other methods we tried. 



C. Summary of known results for the triangular lattice 

In this section we review results on the ground state structure of the supersymmetric 
model on the triangular lattice that were previously established. For further details we refer 
the reader to the original publications. As was mentioned in the introduction a lower bound 
on the number of zero energy states was obtained in [33] by the numerical computation of the 



Witten index (1.2) for the MxN triangular lattice with periodic boundary conditions applied 
along two axes of the lattice. The numerical results immediately indicate an exponential 
growth of the absolute value of the Witten index. To quantify the growth behavior, the 
largest eigenvalue A at of the row-to-row transfer matrix for the Witten index on size MxN 
was determined. For the triangular lattice this gives |33| 

|A| ~ 1.14 , arg(A) ~ 0.18 X TT 

leading to a ground state entropy per site of 

fl > ^ log \Wm,n\ ~ log |A| ~ 0.13 . 

The argument of A indicates that the asymptotic behavior of the index is dominated by 
configurations with filling fraction around u = F/{MN) = 0.18. 

There are two main results for bounds on the total number of zero energy states for the 
triangular lattice. Both results were obtained by considering the (co)homology problem for 



the independence complex associated to the triangular lattice (see section IBl). The first 
result, obtained by Jonsson [28j, proves that zero energy ground states exist in a certain 
range of filling. He constructs a certain type of non-trivial homology element called cross- 
cycles. The size of a cross-cycle refers to the number of occupied sites in the homology 
element. If a cross-cycle of size k exists, it follows that the homology for k particles is non- 
vanishing. Using the relation between (co)homology elements and quantum ground states 
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it follows that there exists at least one zero energy state with k particles. Jonsson obtains a 
bound on the size of the cross-cycles and this results in a bound on the particle numbers for 
which the homology is non- vanishing. That is, there is a set of rational numbers r, such that 
there exist cross-cycles of size rN, where is the number of vertices of the two-dimensional 
lattice. For the triangular lattice it is found that r G [|, |] fl Q. In other words, there are 
zero energy ground states in the entire range between 1/7 and 1/5 filling. 
Let us give the specific form of a cross-cycle z of size k: 

• z = Y\d=i{\'^i) ~ l^i)) such that 2; is a state in the Hilbert space with k particles, that 
is the ai and hi obey the hard-core condition. 

• Furthermore, there is at least one configuration in z such that all sites in the lattice 
are either occupied or adjacent to at least one occupied site. This is called a maximal 
independent set. 

• Finally, is adjacent to h^. 

Note that, in this case, we consider the homology and not the cohomology. It is easily 
verified that z belongs to the kernel of Q\ since Q'^ gives zero on each term in the product: 

Qt(|a,) - = 10) - 10) = 0. 

The latter two conditions ensure that z is not in the image of Q'^: the second condition 
ensures that there is no site c such that Q^c) riiLilki) ~ 1^*)) — ^ ^"^^ ^^e third condition 
ensures that \aj)\hj) Y\i^j{\ai) — \hi)) violates the hard-core condition. 

Clearly, the latter two conditions, combined with the hard-core condition, impose certain 
bounds on the size of a cross-cycle. For the triangular lattice the size of the cross-cycles is 
at most a fifth of all the sites in the lattice and at least a seventh |28] . 

The second result that imposes a bound on the (co)homology on the triangular lattice, 
was obtained by Engstrom [9j. He finds an upper bound to the total dimension of the 
cohomology for general graphs G using discrete Morse theory. He proves that if G is a graph 
and D a subset of its vertex set such that G\D is a forest, then dim('HQ''') < |Ind(G'[D])|. 
Here ^ is the cohomology of Q on the Hilbert space spanned by all hard-core particle 
configurations with n particles on the graph G and Ind(G[D]) is the Hilbert space spanned 
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by all hard-core particle configurations on the subset D. From this theorem it follows that 
finding the minimal set of vertices that should be removed from G to obtain a forest, gives 
an upper bound on the total dimension of "Hg and thus on the total number of zero energy 
ground states for the supersymmetric model on the graph G. For the triangular lattice of 
size 2m x n the upper bound was found to be approximately 0™", with (f) = ^(1 + V^), the 
golden ratio. 



II. RESULTS ON GROUND STATES OF THE M xN TRIANGULAR LATTICE 



In this section we discuss what we can learn about the ground state degeneracy on the 
triangular lattice in general from our numerical results. The total number of zero energy 
states in each particle number sector was computed numerically for lattices with up to 54 
sites (see appendix [B]) . We compare these results to the upper and lower bounds established 
by Engstrom and van Eerten, respectively. We also compare the fillings for which we found 
zero energy states to the bounds obtained by Jonsson. 

Clearly, one has to be careful drawing conclusions for the full two dimensional triangular 
lattice from relatively small systems. We have decided to exclude the ladder geometries of 
size 2 X L and 3 x L in this analysis for the following reasons. First of all, in these geometries 
the particles cannot hop past each other due to the nearest neighbor exclusion. One readily 
checks that the exclusion rule constraints the configurations to a maximum of one particle 
per rung. Second of all, the lower bound obtained from the Witten index for these geometries 
does not agree with the asymptotic value for the triangular lattice to within the error bars 

The results for the ladder geometries of size 4xL, 5xL, 6xL and 7 x L are shown 



in figure II. 1 A fit to the data suggests that the total number of ground states on the 



triangular lattice of size M x N grows as 

iVgs ~ 1.15*^^. (II. 1) 

This result is in agreement with the upper and lower bounds given by (-\/0)*^^ ~ 1.27*^^ 
and 1.14^^^, respectively. 

We now turn to the range of fillings for which zero energy states exist. As explained in 



section IC, Jonsson obtained a minimal range of fillings for which zero energy states exist. 
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Figure II. 1: We plot the total number of ground states for ladder geometries of size 4 x L (circles), 
5 X L (diamonds), 6 x L (squares) and 7 x L (triangles) on a logarithmic scale as a function of the 
total number of sites. The drawn, gray line is a fit to the data. The dotted, blue line is the lower 
bound and the dashed, red line is the upper bound. 

We now present a maximum on this range. 

Lemma. The cohomology of the M x N triangular lattice is trivial for all fillings z/ > 1/4 
and V < 1/8. 

Proof. We divide the M xN triangular lattice into two sublattices, each consisting of a set of 
diconnected chains. More precisely, sublattice 5*1 consists of all sites with coordinates (m, 2n) 
and sublattice 5*2 consists of all sites with coordinates (m, 2n + 1), with m G {0, . . . , M — 1} 
and n G {0, . . . , \_N/2\ — 1}. We now show that the cohomology of Q2 is trivial for all fillings 
outside the range between 1/8 and 1/4. 

Let us first consider the case that none of the Si sites are occupied. It follows that the S2 
sublattice is a collection of periodic chains. For these chains the cohomology is non-trivial 
at 1/3 filling only. This results in an overall filling of 1 particle per 6 sites. It follows that 
non-trivial elements of Tig at other fillings must have particles on the sublattice. 

The highest possible overall filling is found to be 1/4. On the one hand, maximizing 
the density on Si, which is 1/2 filled, results in an overall filling of 1/4, since all sites on 
the 5*2 sublattice must be empty due to nearest neighbor exclusion. On the other hand, an 
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effective cliain of lengtli 2 on tlie 5*2 sublattice lias a ground state witli 1 particle, that is, 
at 1/2 filling. However, an effective open chain on the 5*2 sublattice always has a blocked 
site at each end, leading to an effective filling of 1/4. One readily checks that maximizing 
the overall filling by occupying the Si sublattice such that the 5*2 sublattice is effectively a 
collection of short chains also leads to a maximal filling of 1/4. 

The minimal filling for which Hq^ is non-trivial is obtained by blocking all 5*2 sites by 
occupying a minimal number of sites in 5*1. One quickly shows that this leads to an overall 
filling of 1/8. 

We conclude that "Hqj is non-trivial for all fillings between 1/8 and 1/4. Using the fact 
Tig C l-iQ^iT-iQ^) ^ 'Hq2) establish this as the maximal range in which zero energy states 
can occur. □ 



Numerically we observe zero energy states for fillings v G [1/7, 1/5] fl Q (see figure II. 2). 
We note that the cross-cycle at 1/7 filling obtained by Jonsson requires systems of size 
M X N, such that both M and are multiples of 7. We therefore only observe 1/7 filling 
for the system of size 7x7. Interestingly, zero energy states at 1/5 filling are observed even 
for systems with periodicities that are not compatible with the cross-cycle at 1/5 filling. 
An example is the system of size 4 x 10. Finally, we note that although the data seems to 
suggest that there are no zero energy states outside the range of fillings obtained by Jonsson, 
we have insufficient data to go beyond speculation. 

Finally, we implemented translational invariance to obtain the ground state degeneracy 
per momentum sector. We resolved for momenta along both periodicities, u = (0, M) and 
V = (A^, 0). The eigenvalues of translations along u and v we write as ty = e^'^'^^v/'^ and 
tx = e^'^'^^^l^ ^ respectively. The data clearly reveals a two dimensional flatband dispersion, 
that is, we typically find zero energy states for all momenta {kx-, ky). An example is shown 
in table [T| We observe this property for all system sizes, with the exception of the 3 leg 



ladder of size 3 x L with L odd. We discuss this special case in more detail in section |III B 
We conclude that the triangular lattice of general size M x N will typically exhibit a two 
dimensional flatband dispersion at least for a certain range of fillings. This feature was also 
observed for the square lattice 
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Figure II. 2: We plot the fillings for which zero energy states exist as a function of the total number 
of sites for ladder geometries of size 4 x L (circles), 5 x L (diamonds), 6 x L (squares) and 7 x L 
(triangles). Open symbols indicate the possible fillings for a given system size. Filled symbols 
indicate the fillings for which zero energy states exist. The red and blue lines are lower and upper 
bounds on the fillings, respectively. The dotted lines indicate the range 1/8 < z/ < 1/4, outside this 
range no zero energy states exist. The dash-dotted lines indicate the range 1/7 < i' < 1/5, inside 
this range zero energy states exist for the full triangular lattice. 

III. SUPERFRUSTRATION FOR 3 AND 4 LEG LADDERS 



In this section we focus on the ladder geometries of size 3 x L and 4 x L. For the 4 leg 
ladder we present a slightly sharper upper bound to the total number of ground states and 
we conjecture a sharper bound on the range of filling for which zero energy states exist. We 
then turn to the 3 leg ladder, for which we present a variety of analytic and numerical results. 
Interestingly, the ground state structure of this system is very different for odd and even 
lengths, L. For L even we find the total number of ground states for each particle number 
sector analytically by solving the cohomology problem, for odd L this result is still lacking. 
We do, however, present a rigorous lower bound on the cohomology and conjecture that it 
is exact [32]. Furthermore, numerical computations for the ground state degeneracy in each 
momentum sector clearly show a flatband for L even for momenta in both the vertical and 
the horizontal direction. For L odd we do not observe a flatband. 
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Table I: We show the ground state degeneracy per momentum sector for the triangular lattice 
of size 6x9. There are zero energy states in the sector with 9 (left) and 10 (right) particles. 
The momentum sectors are labeled by kx and ky as defined in the text. It is clear that there 
are zero energy states in all momentum sectors for both particle number sectors, indicating a two 
dimensional flatband dispersion. 
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A. 4 leg ladder 

For the 4 leg ladder we have not been able to solve the full cohomology problem. However, 
(for L even) we can compute the cohomology of T-Lq^, where Q2 acts only on the sites of 
even rungs and Qi acts only on sites of odd rungs, such that Q — Qi + Q2. This gives an 
upper bound on the full cohomology solution, since T-Lq C %Q^{T-iQ^) C I-Lq^. 

Let us define 5*1 and 5*2 as the sublattices consisting of the sites of odd and even rungs, 
respectively. To compute "Hqj we consider all configurations on sublattice Si such that the 
cohomology of Q2 is non-trivial. Note that the sublattice 5*2 is a collection of periodic 4 
site chains. Particles on block sites on these chains effectively cutting them into smaller 
open chains. For a given even rung in 5*2 we thus have four possbilities: 1) the adjacent 
rungs are empty, leading to a 4 site periodic chain, or the configuration on the adjacent S\ 
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Figure III.l: We show an example of a configuration on sublattice Si that leads to the 4 possible 
effective chains on sublattice 82- We denote sites on sublattices Si and S2 by triangles and circles, 
respectively. Filled triangles are occupied Si sites, the open triangles are empty sites. On the S2 
sublattice we indicate sites that must be empty by nearest neighbor exclusion by the red crosses. 
Finally, the numbers labeling the S2 rungs refer to the 4 possibilities described in the text; 1) 4 site 
periodic chain, 2) 2 site open chain, 3) an isolated site, 4) entirely empty. 



rungs is such that the 5*2 rung is effectively 2) a 2 site open chain, 3) an isolated site, 4) 



entirely empty (see figure III.l). In all situations there is one non-trivial element in Hq^, 
except in the third situation; when the S2 rung is effectively an isolated site the cohomology 
vanishes. It thus quickly follows that a complete basis of representatives of is given by 
all configurations on Si such that there is no 5*2 rung that is effectively an isolated site. 

A transfer matrix generating the allowed configurations is easily constructed. We define 
the activity on sublattice Si as e'^' and write y = e^^ and x = e^"^ . The transfer matrix for 
adding two rungs can then be written as 

/ 2 2 \ 

I X xy xy xy xy y y \ 

X Q xy Q y y"^ y"^ 

X y xy y"^ y"^ 

X y xy y"^ y"^ 

X xy y y"^ y"^ 

1 y y y y y^ y^ 
V 1 y y y y J 

4 (— 1, —1))]") correctly reproduces the Witten index results for 



Tl{x,y) 



(III.l) 



We readily verify that Tr([r, 
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Figure III. 2: We plot the total number of zero energy states for the 4 x L ladders on a logarithmic 
scale. The drawn, magenta line is a fit to the data. The dotted, blue line is the lower bound and 
the dashed, red line is the upper bound obtained by Engstrom and the dash-dotted, orange line is 
the upper bound from the dimension of 'Hq2- 

ladders of length 2n. The dimension of 'Hq2 for a ladder of length 2n is given by 

dim(?^Q,(L = 2n)) = Tr([T|(l, 1))]") = J] = 6" + (-2)" + (-1)", (III.2) 

i 

where Aj are the eigenvalues of the transfer matrix, which are found to be 6, —2, —1, 0, 0, 0. 
It follows that for large n, we obtain a new upper bound to the total number of zero energy 
states. For ladders of even length L we find it to be (6^^'^)^'^ ~ 1.25^^, which is slightly 
sharper than the upper bound by Engstrom {^f^^ ~ 1.27^^). The result can be extended 
to ladders of odd length by constructing a transfer matrix for adding two 5*2 rungs and one 
5*1 rung. This is straightforward, but somewhat tedious. We merely state the result 
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5x X y X y x y x y xy xy 

x^ 2xy 2xy xy"^ xy"^ 

2xy 2xy xy'^ xy'^ 

x^ 2xy 2xy xy"^ xy^ 



X 



2xy 2xy xy xy 



X xy xy xy xy y^ y^ 
y X xy xy xy xy y^ y^ j 



(III.3) 



From this we find 

dim(?^Q,(L = 2n + 1)) = Tr(T|(l, 1)[T|(1, 
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7 



-1) 



n-l 



(III.4) 



which leads to the same asymptotic bound on the total number of zero energy states as for 
the even length ladders. 



In figure III.2 we plot the numerically computed total number of zero energy states 
together with the new upper bound and the two known bounds. We see that the numerical 
results are consistent with these bounds. A fit to the numerical data suggests that the total 
number of ground states goes as ~ 1.16"^^. It is clear that the new upper bound still leaves 
much room for improvement. 

The cohomology solution of Q2 contains more information: it gives an upper bound on 
the number of zero energy states at each particle number. One easily checks that the non- 
trivial elements in H-q^ have at least 1 and at most 2 particles per 2 rungs. This proves that 
there are no zero energy states for all fillings z/ > 1/4 and v < 1/8. In other words, we find a 
strict bound on the range of filling for which zero energy states can exist: 1/8 < z/gs < 1/4. 
However, the fact that (XimiJiQ^ ^ suggests that the actual range will be significantly 
smaller. In figure III.3| we show the numerical results for the range of fillings for which zero 
energy states exist. The results clearly obey the strict bound. Apart from the strict bound 
we have also indicated the minimal range, 1/7 < z/ < 1/5, for which zero energy states exist 
as shown by Jonsson. This result, however, is valid provided that the triangular lattice is 
wrapped on a large enough torus. For the 4 leg ladder it is easily seen that the periodicity 
V = (0,4) does not accomodate the cross-cycles with 1/5 filling or 1/7 filling. We do find 



cross-cycles with 1/6 filling and with 3/16 filling (see figure III. 4). By combining their unit 
cells, we can construct cross-cycles for all fillings in the range 1/6 < u < 3/16 provided the 
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Figure III. 3: We plot the fillings for which zero energy states exist as a function of the total 
number of sites for the four leg ladder. Open circles indicate the possible fillings for a given 
system size. Filled circles indicate the fillings for which zero energy states exist. The red and blue 
lines are lower and upper bounds on the fillings, respectively. The dotted lines indicate the range 
1/8 < < 1/4, outside this range no zero energy states exist. The dash-dotted lines indicate the 
range 1/7 < < 1/5, inside this range zero energy states exist for the full triangular lattice. The 
dashed lines indicate the range 1/Q < v < 3/16, inside this range zero energy states exist for the 4 
leg ladder. 

length of the system is long enough. This thus constitutes a strict minimum on the range 
in which zero energy states exist. The numerical data is in good agreement with this bound 
and may even suggest that for large enough L there are no zero energy states outside this 
range. Based on the data we are confident to conjecture that there are no zero energy states 
for fillings z/ < 1/6. For the maximal filling, however, the data is not conclusive and it might 
still be higher than 3/16. 



B. 3 leg ladder 

We now consider the 3 leg ladder. We start with the results that we obtain analytically: 
the Witten index, the solution of the cohomology of Q for L even and various bounds on 
the number of zero energy states and their filling for L odd. 
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Figure III. 4: We show examples of cross-cycles with 1/6 (top) and 3/16 (bottom) filling (see section 



IC). We use the notation introduced in |28]; the black sites form the set {oj} and the white sites 
form the set {&?;}. Both sets form a maximal independent set. The dashed red line serves as a guide 
to the eye, to see that the periodicities of the cross-cycles are compatible with the 4 leg ladder. 
The figure also shows that by concatenating the two patterns we obtain cross-cycles for all rational 
fillings between 1/6 and 3/16 filling provided the length of the ladder is sufficiently long. 



A closed expression for the Witten index of the 3 leg ladder of general length L is readily 
obtained from the transfer matrix that adds a single rung. The transfer matrix reads 



(III.5) 



( 1 z z z\ 
1 2 
1^00 

\ 1 z y 

where we defined the activity z = e^. The Witten index follows from setting the activity to 
-1: 



W{L) = Tr(T3^(-l)) (III.6) 

= Y.^- (111.7) 

i 

= 2^/^ (e"*^^/^ + e'^-^/^) + e"*^-^/^ + e'""^^^ (III.8) 

= 2^/2+1 cos(7rL/2) + 2 cos(7rL/3) (III.9) 

-1)^/22^/2+1 ^ fo^ I g^gj^^ 

-5l for L odd, 
21 



(III.IO) 



where 5^ = 2 if L is divisible by 3 and 5^ = — 1 otherwise (interestingly, this is the same 
5l that Jonsson introduced to write the Witten index for the square lattice [26]). Note the 
starck contrast between ladders of even and odd length. For the former the absolute value of 
the Witten index grows exponentially with system size, while for the latter the lower bound 
on the number of ground states it at most 2. 

We now consider the cohomology problem for L even. We define 5*2 as the subset of all 
sites on the even rungs and 5*1 as the rest. The cohomology of I-Lq^-, where Q2 acts only 
on the 5*2 sites is found from considering all possible configuration on ^i. For empty, S2 
consist of L/2 periodic chains of length 3. Each chain has 2 ground states with 1 particle. 
So we find that the are 2^^'^ non-trivial elements in T-Lq^ with /i = and /2 = L/2. Now 
if we occupy a site on Si it blocks two sites of the 5*2 chains on the two neighboring rungs. 
These rungs are thus effectively isolated sites. The cohomology is then trivial unless these 
isolated sites are also blocked. Note that we can have at most one particle per rung. It 
follows that the cohomology is trivial for all /i > 0, except if /i = L/2. Upon inspection 
one can show that the number of ways to block all S2 sites is: 



2^/2 + ^-^(_l)i/2 



(III.ll) 



Consequently, (III.ll) gives the number of non-trivial elements in I-Lq^ with /i = L/2 and 
/2 = 0. This solves "Hqj- easy to show that T-iQ^i/HQ^) has the same non-trivial elements 
by checking that within I-Lq^ all elements are in the kernel of Qi and therefore not in the 
image of Qi. Finally, since all elements have / = L/2, we immediately find "Hg = "Hq^ (/Hq^)- 
So the number of zero energy states is 



Ngs{L) = 2^/2+1 + ^^(_i)V2 for I gygn^ (III.12) 

All ground states have L/2 = MN/6 particles, where MN = 3L is the total number of 
sites. This corresponds to 1/6 filling. We find perfect agreement with these results and the 
numerical results for ladders of even length up to L = 16. 

For ladders of odd length L we also study the cohomology problem. We make the same 
subdivision in lattices 5*1 and 5*2, but now Si consists of (L — 3)/2 disconnected rungs and 
one pair of adjacent rungs. Following the same reasoning as for L even, we now find that 
Hq^ consists of 2^^^-'^)/^ non-trivial elements with /i = and /2 = (L — l)/2 and 2^^"-'^)/^ -f 
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5l{—1)^^~^'^^'^ non-trivial elements with fi = {L + l)/2 and /2 = 0. Again one readily checks 
that 'Hq-^{'Hq^) has the same non-trivial elements. However, to obtain I-Lq from 'Hq^{'Hq^) is 
now a non-trivial step and turns out to be very challenging. We will comment on this below. 
We first mention a few results that we obtain directly from 'Hqi{'Hq^). Since 'Hq^{'Hq^) 
contains the cohomology of Q, it gives an upper bound on the total number of zero energy 
states: 

iVgs(^) < <lim{HQ,{HQ,)) = 2(^+i)/2 + for L odd. (III.13) 

Furthermore, we find that the cohomology of Q is non-trivial only in the sectors with (L ± 
l)/2 particles. For large L this converges to 1/6 filling. 

In figure [III. 5 we plot the numerical results for total number of zero energy states for 



L odd. We also plot the upper bound given by the dimension of 'Hq-^{'Hq^) and the lower 
bound given by the Witten index. There is a clear substructure with period 3, both in the 
numerical data as well as in the analytic results. We thus extract the asymptotic behaviour 
of the number of ground states by fitting to the running average of the data over three 
consecutive ladder lengths. We find Ngs ~ 1.08^*^. Note that, in contrast, for the 3 leg 
ladder with an even number of rungs we found A^gs ~ 1.12^^"^. This result clearly indicates 
that the last step in the cohomology computation, i.e. finding "Hq from 'Hqi{T-Lq2) is a 
non-trivial step. Finally, we mention that all zero energy states are found in the sectors 
with (L ± l)/2 particles in agreement with the result for 'Hq-^^I'Hq^) ■ We also observe that 
the number of zero energy states with (L + l)/2 particles equals the number of zero energy 
states with {L — l)/2 particles up to ±6l in agreement with the Witten index. 

The numerical results for the number of zero energy states, not only resolve for particle 
number, but also for momentum along the legs and the rungs of the ladder. For the 3 leg 
ladder this data reveals again a very different structure for ladders with even length and 
ladders with odd length. For the ladders with L even the number of zero energy states 
is distributed evenly over all momentum sectors. This implies that there is a completely 
flat dispersion relation in all directions. This flatband property was also observed for the 
supersymmetric model on the square lattice [20] . For the square lattice this could be un- 
derstood from a mapping of the cohomology elements onto tilings |20| [2ll [26]. The unit cell 
of a tiling directly translates into the eigenvalues under translations of the corresponding 
ground state. For the 3 leg ladder we also flnd a correspondence between representatives 
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Figure III. 5: We plot the total number of zero energy states for the 3 x L ladders with odd L on a 
logarithmic scale. The drawn, magenta line is a fit to the running average of the data over three 
consecutive ladder lengths. The dotted, blue line is a lower bound given by the absolute value of 
the Witten index and the dashed, red line is an upper bound on the number of zero energy states 
given by the dimension of 'Hq^{'Hq2)- 



of the cohomology and tilings. An example is shown in figure |III.7 There are two tiles 



and each tile can be neighbored by both tiles (this corresponds to the fact that the ground 
state degeneracy goes as 2^/^). It is easy to see that there are many possible tilings with 
a unit cell equal to the size of the system. Each such tiling thus corresponds to 3L ground 
states, one in each momentum sector. This explains the observed fiatband. Tilings with a 
smaller unit cell correspond to ground states that occur only in certain momentum sectors. 
We have checked that for finite systems a careful analysis of all tilings and their properties 
under translations indeed reproduces the distribution of ground states over momentum sec- 
tors as found in the numerics. We note that just as for the square lattice the correspondence 
between cohomology elements and tilings is not exactly 1-to-l. In fact we find 

N,,{L) = iVtnings(^) + hi-l)"-^', (111.14) 

which is remarkably similar to the result obtained for the square lattice. 

For the 3 leg ladder with L odd we observe the following. When L is a multiple of 3 there 
are zero energy states in the sectors with = 2n — Py = 2nk/3, with A; = 0, 1, 2. When L 
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Figure III. 6: Tilings on the 3 leg ladder. Above we show the 2 allowed tiles in red and green. Note 
that it is the same tile with 2 possible orientations, the third orientation is not compatible with 
the boundary conditions. Below we show an example of a tiling with periodicities u = (0, 3) and 
v = (12,0). 

is not a multiple of 3 all zero energy states have Px = Py = 0. The substructure with period 
3 in the systemsize is again apparent. 

The results presented in this section clearly indicate that the cohomology problem for 
3 X L with L odd is highly non-trivial. Nevertheless, we have been able to make substantial 
progress beyond the results mention above. We believe a complete solution is possible and 
hope to present it in forthcoming work [32j. At this point we present a rigorous lower bound 
on the number of ground states. 

Lemma. The number of zero energy states, Ngs, for the 3 x L triangular lattice with L odd 
obeys: 

ATgs > const. X [(0)^/^]^^ ~ const, x 1.08^^, (III.15) 
with (f) = ^(1 + the golden ratio. 

This lower bound agrees with the exponential growth observed in the numerical data 
and, in fact, we believe that (0)^^^ is the exact base for the exponential growth. We limit 
ourselves to a short sketch of the proof of this lower bound and defer the details to |32] . 

The key ingredient is to compute 'Hq^{'Hq2) for a different choice of the sublattices 5*1 
and 5*2. Instead of the choice discussed above, we now take Si to be one rung and 5*2 the 
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rest of the system. 5*2 is thus an open 3 leg ladder of even length. The cohomology of Q2 
is still easily obtained for this choice of sublattice, however, we pay extra care to choose 
the basis in such a way that the action of Qi is relatively simple. As before we find that 
the number of non-trivial elements in T-Lq^ is 2'^^^^''^'^ + l)'-'^^^^/^, however, we now find 
that 2(^-i)/2 of these elements have /i = and /a = (L - l)/2 and 2(^-i)/2 + 5^(_i){i-i)/2 
of these elements have /i = 1 and f2 = {L — l)/2. As a consequence, the next step in the 
cohomology computation, that is, finding 'Hq^{T-Lq^), is now non-trivial. The upshot is that 
for this construction we immediately have "Hg = 'Hq^{'Hq^). It is possible to express the 
action of Qi on the non-trivial elements of T-Lq^ in a general form, that is for general length 
L. From these expressions, we can identify certain linear combinations of the basis elements 
of I-Lq^ that are annihilated by Qi up to terms that are in the image of Qi. These elements 
are thus in 'Hq^{'Hq^). For L = mod 3 we find that these elements can be written as 
sequences of tiles on 5*2. More precisely, we can identify an element of 'Hqi{'Hq2) with each 
tiling of 5*2 up to cyclic shifts. The rigorous lower bound then follows from counting the 
number of these tilings (note that the reduction due to cyclic shifts is at the worst a factor 
of 2/(L — 1), which does not affect the exponential behavior). The tilings consist of 4 types 



of tiles depicted in figure III. 7 



The counting of the tilings is easily done. We define t{l) as the number of tilings at 
L = 2Z + 1. The counting of t{l) leads to a recursion 



t{l) = t{l - 2) + t{l - 4) + 2t{l - 3). (III.16) 

The characteristic equation has a factor (x^ — x — 1) and the largest eigenvalue is the golden 
ratio = (H-v^)/2. The number of ground states growth as (0^/^)^-^ ~ l.OS^'^ in agreement 
with the numerical results. The resulting value for the ground state entropy per site is 

Sgs/{5L) > log(0)/6 ~ 0.080 . (III.17) 

We conclude that other than for L even, the system hesitates between two values for the 
number of particles in the ground state, L/6 ± 1/2. This gives a cancellation that reduces 
the number of ground states we see for L even, 2^/^, to the smaller number given above. 
Other than is suggested by the Witten index, the growth is still exponential in L. This is 
the second example, after the analysis for the square lattice j2l], where a non-trivial ground 
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Figure III. 7: We show the four types of tiles that are used to construct cohomology elements for 
the 3 leg ladder of odd length, L. Note that the tiles are composed of the two tiles introduced for 
the 3 leg ladder of even length (see figure III.6). 



state degeneracy is derived from a tiling argument. At this instance the analysis is less 
complete though, since the full solution of the cohomology problem is still open. There are 
two observations that make us confident to conjecture that (0)^^^ is the exact base for the 
exponential growth. The first is that it nicely agrees with the fit to the numerical data. 
The second is that the number of tilings for L = 15 and L = 21 seems to reproduce the 
number of zero energy states in the two non-trivial momentum sectors, (27r/3, 47r/3) and 
(47r/3,27r/3), namely t(7) = 10 = 5 + 5 for L = 15 and t(10) = 44 = 22 + 22 for L = 21 (see 
appendix B 1\ . 



IV. COMPLETE SOLUTION FOR THE 2 LEG LADDER 

In this section we discuss the simplest ladder obtained from the triangular lattice by 
imposing the periodicities u = (0,2) and v = {L,0). The resulting ladder geometry is 



depicted in figure IV. 1 This lattice is special because it is invariant under the exchange 
of the 2 sites on a rung. A state on the lattice will have either even or odd parity under 
this transformation. Consequently, this local Z2 symmetry distinguishes an exponential 
number of sectors. This observation will allow us to understand the exponential ground 
state degeneracy. Furthermore, we find a mapping from the ladder to the chain such that 
the entire spectrum can be understood. Remarkably, in certain sectors the spectrum turns 
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Figure IV. 1: The 2 leg ladder of length L, with sites labeled 1, . . . , 2L. 

out to be gapped, in others it is gapless and in yet other sectors we find phase separation. The 
continuum hmit of each of the gapless sectors is shown to be described by a superconformal 
field theory with central charge c—1. 



A. Mapping the ladder to an infinite number of chains 

In this section we present a mapping that maps the supersymmetric model on the ladder 
onto the supersymmetric model of a collection of chains. As a first step we write the three 
possible states on a rung i, that connects the sites 2i — l and 2i, as follows: the empty rung, 
which we denote by |0)j and two states with one particle on the rung: 

l±)i = ^(4-i±4)|o).. (iv.i) 

The next thing to note is that the state is locahzed and carries zero energy. The 
locahzation is due to the fact that there is negative interference for the particle to hop off 
the rung. The state has zero energy because the hopping within the rung precisely cancels 
the potential energy of the state. The effect of these localized zero energy states on the rungs 
is that they effectively cut the ladder into smaller ladders with open boundary conditions. 

Another way to understand this breaking up of the Hilbert space is via the parity operator, 
Tj, that interchanges the sites 2i— 1 and 2i. One easily checks that the Hamiltonian commutes 
with this operator. It follows that the parity on each rung is a good quantum number that 
is conserved under the action of the Hamiltonian. Clearly, the state |— )j has odd parity 
under this tranformation, whereas the states |0)i and have even parity. The Hilbert 
space breaks up into disconnected sectors characterized by the parity of each rung. 

One thing that is important to note is that the rungs neighboring a rung with odd parity 
must be empty due to the nearest neighbor exclusion and must thus have even parity. 

Suppose that the parity of two rungs that are 2 + L_(. rungs apart is odd and that all 
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rungs in between these two rungs have even parity. It follows that there are L+ consecutive 
rungs that can be in either of the two even parity states. We will now show that the Hamil- 
tonian acting on these consecutive rungs with even parity maps to the supersymmetric 
Hamiltonian on the chain of length L+ up to an overall factor of 2. In this mapping the 
empty rung will map to the empty site on the chain and the state will map to occupied 
site on the chain. It is easy to see that the Hilbert spaces are indeed identical. On the chain 
there are two states per site, on the ladder there are two states per rung, furthermore, the 
nearest neighbor exclusion ensures on the ladder that two neighboring rungs cannot both 
be in the |+) state and on the chain that two neighboring sites cannot both be occupied. 

Let us introduce the creation operator c| = {c\^_i + c\^)/V2, which creates the from 
the empty rung. Furthermore, we introduce the number operator -hi = n2i-i + n2i and 
the projection operator P(^j^ = (1 — ri,(j_i))(l — ^(j+i)). Note that the latter is indeed a 
projection operator, because due to the nearest-neighbor exclusion only takes the values 
or 1. When we label the L+ even parity rungs 1 through L+, the Hamiltonian on that 
part of the system reads 

2L+ 2L+ 

i=0 j6(i> i=0 

Expanding these terms and carefully doing the algebra, we can rewrite this Hamiltonian 
in terms of the operators we just introduced which are defined on a chain of sites (we 
deferred the details to the appendix lAl) . What we find is 



H+ = J2Y1 24>£JqAo + E 2Ai> = 2i^chai„, (IV.3) 

j=l i=j±l j=l 

which is precisely twice the Hamiltonian on a chain of sites with an additional empty 
site on each end of the chain. This is equivalent to saying that the chain has open boundary 
conditions. Finally, we note that the Hamiltonian in the sector where all rungs of the ladder 
have even parity maps to the Hamiltonian on the chain with periodic boundary conditions, 
again up to a factor of 2. It follows that in the sector with odd parity on the rungs {li, . . . ,ls} 
and even parity on all other rungs, the Hamiltonian on the 2 leg ladder takes the form 

S Zm + 1 — 1 

^ = 2 E [ E ( E + PU))] ^ (IV.4) 

m=l j=lm. + l «=J±1 
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where we defined Ig+i = h- 

This mapping is very powerful because we have very good understanding of the super- 
symmetric model on the chain [T51ITU] . In the next sections we will first see how we can use 
it to get a closed expression for the number of zero energy states and then how it allows us 
to fully understand the low lying spectrum. 

B. Ground state partition sum 

In the previous paragraph we have seen that the state with odd parity on the rung is a 
dark state: it is localized and decoupled from the rest of the system. In addition it carries 
zero energy. The effect of these zero energy states is that they effectively cut the ladder 
into smaller ladders. These smaller ladders can be mapped to the chain, for which we know 
the number of zero energy states. Remember that the chain of length 1 mod 3 has no zero 
energy states. It follows that in a zero energy state on the ladder the number of rungs 
separating two odd parity rungs cannot be 3/ rungs apart. This is because the even parity 
rungs will map to an open chain of length 3/ — 2 for which there is no zero energy state. 

Carefully counting the number of zero energy states this way, one obtains the ground 
state partition sum as a function of the fugacity z: 

Zl{z) = TrQgz^ 

= z{Zl_2{z) + 2Zl_,{z)) - Al, (IV.5) 

with 

f for L = 3/ and L = 3/ - 1 

Al = < (IV.6) 

I z'+i + z^ for L = 3/ + 1 

Since the number of sites, NM, of the ladder is 2L, it follows immediately from the recursion 
relation that there will be zero energy states for f /{NM) G [1/6, 1/4]. 

To obtain the number of zero energy ground states we solve the recursion relation for 
Zi{z) and set 2; = 1. Using Z2 = 3z, Z^ = 5z and Z4 = Iz^ + we find 

Zl=ziUz = 1) = (-l)^Af/' + {-\)\\2fl^ + \T - (-1)", (IV.7) 

with a = —1, 0, 1 and Ai^2 and A3 the complex and real solutions to —4 + A — 2A^ + A^ = 0, 
respectively. Note that A2 = Ai. For large L, we have 

Zl{z = 1) ^ Ag/^ ^ 1.5214^. (IV.8) 
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It follows that the ground state entropy per site is 

Sgs ^^IuZl^ lnl|214 ^ ^ ^ggg . . . (IV.9) 
Equivalently, the Witten index is obtained by setting z = —1. We find 

H.,= (Zl±!^). + (Zi^,. (iv.lO) 
For large L the Witten index grows as ~ 1.14-^. 



C. Spectrum 

1. Gapped at 1/4 filling 

For L even there are two zero energy states with L/2 particles. The two ground states 
are product states of resonating single particle states (|— )) on every other rung. One of the 
states has particles resonating on all odd rungs, the other on all even rungs. For / even 
(odd) translation by two rungs maps each ground state to — (+) itself, so the eigenvalue of 
the translation operator obeys: — —1 ior f even and — 1 for / odd. It follows that the 
translationally invariant ground states have momenta: 

Pi — 7r/2 and p2 — 37r/2 for / even 
Pi = and P2 = TT for / even. 

Note that quarter filling is the highest possible density for this system due to the nearest- 
neighbor exclusion. 

For L odd the maximal number of particles in the system is (L — 1)/2. We find that there 
are 2L zero energy states with / = (L — l)/2. Pictorially, we can write these states as 

/ /-I 

i^)=io)nio-) i5)=(i+o)-io+))/^nio-), (iv.ii) 

fc=i fe=i 

where denotes an empty rung and ± denote the single particle states with even and odd 
parity, respectively. Clearly, there are L states of each type. These L states are related by 
translations by one rung. It follows that the translationally invariant ground states have 
momenta pk — 2'Kk/L with k — 0, . . . L — 1 and there are 2 in each momentum sector. The 
two types of states occur in different parity sectors: the states of type A occur in the sector 
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with / rungs with odd parity and the states of type B occur in the sector with / — 1 rungs 
with odd parity. Note that in the latter sector, there can be 4 consecutive rungs with even 
parity, which can be mapped to a 2-site open chain. Indeed the ground state on a 2-site 
open chain is (| • o) — | o •))/^/2. 

At quarter filhng there is a gap to charge-neutral excitations. The excitations at quarter 
filling follow from changing one of the single particle states from |— ) to 1+). This costs an 
energy A = 2. For L even there are L such excitations, whereas for L odd there are — 2L 
such excitations^. Consequently, the spectrum is gapped to chargeless excitations, the gap 
is A = 2 and the first band above the ground state is completely flat. Two 1+) -excitations 
obey a certain exclusion statistics, so the distribution of states over the momentum sectors 
is not uniform, but the band is still fiat &t E = 2A = 4. Continuing this way, the highest 
energy state is reached aX, E — L ioi L even (for L odd the highest energy state is found to 
have E^ L + 1). 



2. Gapless at 1/6 filling 

At 1/6 filling the number of ground states is exponential in the length. In the sector with 
even parity on all rungs there are two zero energy states with f = j for L = 3j corresponding 
to the two ground states of the periodic chain of length L = 3j. Similarly, there is one zero 
energy states with f — j ior L — 3j ± 1 corresponding to the ground state of the periodic 
chain of length L — 3j ± 1. In the sector with j rungs with odd parity, the zero energy 
states can be depicted as 



|gs) = < 



- 0) nU io - 


-0) 


for L — 


3j, 


00-) m=i io - 


0) 


for L = 


3j + 1, 


-o)ni=ilo- 


0) 


for L = 


3j - 1. 



(IV.12) 



These states are clearly product states, which implies a finite correlation length and thus 
a gap to charge-neutral excitations. In contrast the ground states in the sector with even 

parity on all rungs, correspond the the ground states on a chain, which are known to be 
critical. It follows that the spectrum is gapless in that sector. In fact, the supersymmetric 

^ Note that for L odd there is the extra possibihty of leaving the |— ) states intact, but exciting the length 
2 ladder to its first excited state, however this costs an energy A = 4. 
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model on the chain is well understood and the low-energy spectrum is known to be described 
by an A/" = (2, 2) superconformal field theory with central charge, c = 1 [T^ [TU] . In (TH] a 
complete identification between states in the spectrum of chains of finite length and fields 
in the continuum theory is given. This identification allows us to compute the gap scaling 
for the 2 leg ladder in each parity sector. Interestingly, we will find that the lowest non-zero 
energy state is not in the sector with even parity on all rungs. 

The finite-size scaling of the energy of a gapless system depends on the boundary condi- 
tions: 

E,hain = bnEscFTVp/L + Oil/L^), (IV.13) 

with 6 = 1 for open and 6 = 2 for periodic boundary conditions. For the supersymmetric 
model on the chain the Fermi velocity was found to be vp = 9^/3/4: [19]. For the periodic 
chain the first excited states have -Escft = 2/3 for L = 3j and -Escft = 1/3 for L = 3j ± 1. 
Using the fact that H = 2ifchain5 we find that the finite-size energies of these states on the 
ladder are 

^ , f 6nV^/L for L = 3j ^ ^ 

^ladder = TtEsCFT^VS/L ={ ' {IV .U) 

I 37iV3/L for L = 3j± 1. 

For the open chain with L = 3j and L = 3j — 1 we find that -Escft = 1 for the first excited 
state. However, for the open chain of length L = 3j + l there is no zero energy state. Instead 
the lowest energy state at / = j has -Escft = 1/3 and this state has a superpartner with 
the same energy at / = j + 1. Using the above, we find that the finite-size energies of these 
states on the ladder are 

7rEscFT9y3 I 97rv^/(2L) for L = 3j and L = 3j - 1 

-^ladder = ^7 = < ^ (IV.15) 

2^ 1^ 3nV3/{2L) for L = 3j + 1. 

Comparing all the above energies, we find that the lowest non-zero energy state of a ladder 
of length L corresponds to a state on an open chain of length 3m + 1 = L — x, with x = 3, 5, 7 
for L = 3j + l,3j,3j — 1, respectively. The scaling of the energy of the first excited state 
thus reads 



37rv^ 



3nV3 
2(L-3) 



for L = 3 j + 1 



2(3m + 1) 



for L = 3j (IV. 16) 

for L = 3j — 1 



2(L-5) 
I 2(L-7) 
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The corresponding states are found in the sectors with odd parity on the rungs {/}, {/, / + 2} 
and {1,1 + 2,1 + 4}, respectively, where / G (1, . . . , L). Note that the states on the chains 
have / = j — 1, j — 2, j — 3, respectively, so that the total number of fermions is always 
f = j- There are no zero energy states in these sectors. Finally, it is clear that the first 
excited state is L-fold degenerate: one for each momentum sector pk = 2TTk/L. 

Numerical results for ladders up to L = 20 confirm these observations. In particular, a 
finite-size scaling analysis shows that the energy of the first excited state is nicely fitted with 
ai/L + a2/L^ + a^/L^. We extract Oi ~ 8.17 for L = 3j + 1 and Oi ~ 8.15 for L = 3j, which 
is good agreement with the theoretical value 37r-\/3/2 ~ 8.162. The result for L = 3j — 1 



(ai ~ 8.42) is less accurate. Indeed, from (IV.16) it is clear that for this length the finite-size 
corrections are the strongest. 



3. Phase separation at intermediate fillings 

We can continue the logic of the previous section to infer the energy of the first excited 
state for fillings // (2L) G [1/6, 1/4]. Starting from the 1/6 filling side, there are two obvious 
possibilities for the first excited state at higher densities: the density is increased either by 
having more particles in the part of the system where all rungs have even parity, or by going 
to a parity sector where for a larger part of the ladder even and odd parity rungs alternate. 
From the continuum theory we know that in the first case the energy scales parabolically 
with the fermion number: 

£^scft(/) = IP - If, (IV.17) 

where / = /chain ~ /gs ~ 1/3. For the ladder we thus have that the energy of such states with 
f = j + j' particles scales as 



3nV3 (\ 9 



L — X 



Q + (IV.19) 



where x = 3, 5, 7 for L = 3j + 1, 3j, 3j — 1, respectively. 

In the second case, one can show that the first excited state with f = j + j' particles 
corresponds to the lowest energy state of a chain of length 3m' + 1 = L — 6j' — x with m' 
fermions for j' even, and x as before. For j' odd we have 3m' + 1 = L — 6(j' — 1) — x with 
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m' + 1 fermions. It follows that the energy of this state is 



for j' even and similarly for j' odd with j' replaced by j' — 1 everywhere. It is readily checked 



,-/2 

"^2 



that for j' > 0, this energy is always smaller than the energy for the first case (IV. 18). 

It follows that at some intermediate filling between 1/6 and 1/4, the lowest energy states 
show phase separation, where part of the ladder is in the |0 — ■ ■ ■ — 0) phase and the rest is 
in the 1/6 filled critical phase. So starting from 1/6 filling, both charged and charge-neutral 
excitations above the ground state manifold are gapless. We know, however, that at 1/4 
filling, there is a gap to charge-neutral excitations. The question is thus at what filling this 
gap opens up. Suppose the density is 1/4 — e in the thermodynamic limit. We know that 
the first excited state shows phase separation. The energy of the first excited state will scale 
as l/l, where / is the length of the part of the system that is in the 1/6 filling phase. One 
easily checks that / oc eL. It follows that for any finite e the gap collapses. 

Finally, we note that using parity the gapless phase separates into an exponential num- 
ber of sectors^. In the continuum limit, each of these sectors is described by a (sum of) 
superconformal field theories with central charge c = 1. 



D. Relation to other superfrustrated ladders 

In this section we briefly discuss two other ladder geometries that share certain properties 
with the 2 leg ladder obtained from the triangular lattice. 



The first is the octagon-square ladder depicted in figure |IV.2[ which can be obtained 
from the octagon- square lattice by imposing doubly periodic boundary conditions. The 
solution to the cohomology problem for the supersymmetric model on the octagon-square 



^ To obtain an estimate on the number of gapless sectors in the continuum we proceed as follows. As a 
sufficient condition for the sector to be gapless, we require that there is at least one part of the system, 
say Lg consecutive rungs, with even parity, such that Lg goes to infinity as L goes to infinity. To make 
things concrete, let us take Lg = L/2. The parity of the rungs in the rest of the system then merely has 
to be such that a zero energy state exists. From the ground state partition function, |IV.8| we easily see 
that the number of sectors that obey these conditions can be estimated as ~ l.S'^^^' = 1.5^/"^. This 
shows that there is indeed an exponential number of gapless sectors. 
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^^^^^^^^ 

Figure IV. 2: We show the octagon-square ladder (above) and the zig-zag ladder (below). In the 
octagon-square ladder the vertical links are referred to in the text as links across the plaquette. 
This example shows a ladder of 4 plaquettes. In both cases periocid boundary conditions in the 
horizontal direction are implied. In the zig-zag ladder the red line serves as a guide to the eye to 
see that the zig-zag ladder is a chain with nearest and next-nearest neighbor interactions. 

lattice is known for general boundary conditions [14] . For the octagon-square ladder with 
plaquettes, the number of zero energy states is 2^ -|- 1. This exponential ground state 
degeneracy can be understood in terms of a local symmetry, just as for the triangular ladder. 
Let us label the sites around a plaquette in a clock wise direction starting at the leftmost site 
as i,i + + 2, i + 3. It follows that the sites i + 1 and i + 3 are connected by the link across 
the plaquette. It is clear that the ladder is invariant under the exchange of these two sites. 
Moreover, the odd parity single particle state, |— )j = (c|^i — c|^3)/-\/2|0), is again a dark 
state with zero energy. So having odd parity on a plaquette effectively cuts the ladder into 
an open ladder. Finally, there is again a mapping to the chain Hamiltonian if one restricts to 
the sector with even parity on the plaquettes [lOj. There are three main differences. First, 
there are zero energy states for all parity sectors. This can be seen as follows. The odd 
parity plaquettes cut the ladder into smaller ladders of lengths Ij = mod 3 for all j. It 
follows that there is one zero energy state in all parity sectors, except in the sector with 
even parity, where we effectively have a periodic chain of length 3A'^ which has 2 zero energy 
states. Note that there are 2^ parity sectors and thus 2^ -|- 1 zero energy states. Second, 
all zero energy states have filling 1/4. Remember that the ground states on chains of length 
/ = mod 3 have 1 particle on every 3 sites, which corresponds to 1 particle per plaquette. 
Clearly, the plaquette with odd parity also have 1 particle per plaquette. Therefore, all 
ground states have 1/4 filling. Third, and most importantly, the spectrum is not described 
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by a conformal field theory. This can be seen by carefully carrying out the mapping to the 
chain. One then finds that the effective Hamiltonian is the Hamiltonian for the staggered 
chain, H = {Q\, Q\}, with Qx = J2j staggering Xj = \pl for j = mod 3 and 

\j = 1 otherwise. This model was studied in detail in Refs. [TTl [T2l [2T] and has a gapped 
spectrum, unless the boundary conditions are such that they allow for a massless kink. A 
massless kink exists on open chains of length / = mod 3 and has a 1/P dispersion. It 
follows that certain sectors of the octagon-square ladder are gapless, but not conformal, and 
others are gapped. 

The second model to which our 2 leg ladder bears certain similarities is the zig-zag ladder 



(see figure IV.2). This ladder is obtained from the square lattice by imposing periodicities 
u = (1,2) and v = {L,0). The similarities are not as striking as those observed for the 
octagon-square ladder, but we believe that they can be quite useful in trying to get a better 
understanding of the zig-zag ladder. We find that there are two important similarities. 
First, there are similarities in the ground state degeneracy. The cohomology problem for 
the square lattice with doubly periodic boundary conditions was solved for a large class of 
periodicities [20ll2^ . It was found that for the zig-zag ladder there is an exponential number 
of zero energy states that occur in the range of fillings between 1/5 and 1/4. The counting 
of the ground states is rather non-trivial and can be formulated as a tiling problem. Second, 
there are similarities in the spectrum. The spectrum of the zig-zag ladder was investigated 
and it was found that there are gapped as well as gapless sectors [TE\ [20] • In particular, 
the spectrum was found to be gapped at 1/4 filling, and the 4 zero energy states are simple 
product states with a particle resonating on every fourth diagonal rung. This is very similar 
to the 1/4 filled ground states of the 2 leg triangular ladder (although in this case there are 
only 2). For the zig-zag ladder is was found that the spectrum is also gapped at 1/5 filling. 
The gapless phase was observed at intermediate fillings. In particular, 2/9 filling seems to 
play an important role. We think that the ladder considered here can shed light on the 
physics of the zig-zag ladder. What we observe here at 1/6 filling may be a cartoon version 
of what is going on in the zig-zag ladder at 2/9 filling. A promising route could be to try 
to perturb the ladder considered here, in such a way that the dark states become mobile. 
Another possibility is to study the zig-zag ladder in the limit of zero staggering on every 
third site of the zig-zag chain. In this limit the spectrum bears very strong resemblances to 
that of the ladder studied here. In fact, it also falls apart into many sectors, some gapped. 
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some gapless, with the gapless sectors described by (sums of) superconformal field theories 
with central charge c = 1. 

V. CONCLUSIONS 

In this paper we combined numerical and analytical techniques to investigate the ground 
state structure of the supersymmetric model on the triangular lattice with doubly periodic 
boundary conditions. Previous studies showed that this system has an exponential ground 
state degeneracy, leading to an extensive ground state entropy jOl |33]. Furthermore, zero 
energy states were proven to exist between 1/7 and 1/5 filling |28| . Our numerical studies 
confirm these results. The two main results we obtain in this work for the full 2D triangular 
lattice are an analytic upper bound on the range of fillings for which zero energy states can 
exist, namely between 1/8 and 1/4 filling, and the numerical observation that the system 
typically exhibits a two dimensional fiatband dispersion, which we believe to persist for the 
full 2D model. 

We also studied 3 ladder geometries in more detail and presented a variety of new results. 
For the 4 leg ladder we present a slightly sharper upper bound to the total number of 
ground states and we conjecture a sharper bound on the range of fillings for which zero 
energy states exist. For the 3 leg ladder we find that the ground state structure is very 
different for odd and even lengths, L. For L even we find the total number of ground 
states for each particle number sector analytically by solving the cohomology problem, for 
odd L this result is still lacking. We do, however, present a rigorous lower bound on the 
cohomology, which is in perfect agreement with the numerical data, and conjecture that it 
is exact [32j. Furthermore, numerical computations for the ground state degeneracy in each 
momentum sector clearly show a fiatband for L even for momenta in both the vertical and 
the horizontal direction. For L odd we do not observe a fiatband, instead the ground states 
all have = 2it — py = 27[k/3, with k = 0,1, 2. Interestingly, many of the results for the 
3 leg ladder can be understood from a mapping between ground states and tilings. This is 
thus another example where a non-trivial relation between tilings and ground states of the 
supersymmetric model is observed fIU\ [211 [2S1 

Finally, we discuss the full solution for the 2 leg ladder. For this system the ground state 
degeneracy can be understood in terms of a local Z2 symmetry of the lattice. Furthermore, 
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we find a mapping from tlie ladder to tlie cliain sucli tliat tlie entire spectrum can be 
understood. Remarkably, in certain sectors the spectrum turns out to be gapped, in others 
it is gapless and in yet other sectors we find phase separation. The continuum limit of each 
of the gapless sectors is shown to be described by a superconformal field theory with central 
charge c = 1. We argue that this model can shed light on the rich physics observed for the 
supersymmetric model on the zig-zag ladder [T8l |2Q] . 
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Appendix A: Mapping of 2 leg ladder Hamiltonian to the chain Hamiltonian 

In this appendix we give some details of the computation that leads to the identification 
of the Hamiltonian on the 2 leg ladder in the even parity sector with the Hamiltonian on 
the chain. Let us introduce the creation operator cj = (c2j_i + C2j) / -\/2, which creates the 
from the empty rung. Furthermore, we introduce the number operator fii — n2i-\ +n2i 
and the projection operator P^j) = (1 — 7T,(j_i))(l — 77,(^+1)). Note that the latter is indeed a 
projection operator, because due to the nearest-neighbor exclusion rii only takes the values 
or 1. When we label the L_|_ even parity rungs 1 through L_|_, the Hamiltonian on that 
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part of the system reads 

2L+ 2L+ 

Expanding these terms, we can rewrite this Hamiltonian in terms of the operators we just 
introduced which are defined on a chain of L+ sites. We first rewrite the latter term. 

Hpot = Y P(2j-1) + P(2j) 
L+ 

= Yy^^ ~ + (1 - ^2j-l)](l - ^2(j-l)-l - ri20-l) + n2(j-l)_in20_l)) 

x(l - n2(j+i)-i - n2{j+i) + n2(j+i)-in2(j+i)) 

L+ 

= J^(2 - n2j - n2j_i)(l - ^20-1)-! - ^^2(i-i))(l - n2(j+i)-i - n2(j+i)) 

= E(24>-^.), (A.2) 

i=i 

where we have used the fact that njUj^^Qy-^ to j — ^ hardcore fermions in the third 

and the last fine. To rewrite the kinetic term we first note that 

P{2j)clj = (1 - n2j-i)P{j)clj and 

C2jP{2j) = C2j{l - n2j-l)P{j) = C2jP{j), 

where, in the last step, we use the fact that for hardcore fermions C2jn2j-i — 0. It follows 
that we can write the kinetic term as 

L+ 

Hkin = y^[-P(2j)4,C2j-l-P(2j-l) + P(2j-1)4,-1%-P(2j) 

+ (^(2j)4j + P{2j-l)4j_i){c2i-lP(2i-l) + C2iP{2i))] 

i=j±l 

= J^IAj)^^ ~ ^2i-i)4i%-iAi> + Ai>(i - ^2i)4i-iC2iAi> 

+ XI (^0)(l - ^2j-l)4j + Aj>(l - ^2j)4j-l)(c2i-lAi> + C2i^(i>)] 

^+ 

= E[A.)(4,%-i + 4,-iC2i)Po> + E A.>^^c]x/2qP(,)], (A.3) 

j=l i=J±l 
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where in the last hne we used the fact that n2j-icljC2j-i — cljn2j-iC2j-i = in the first 
term. In the second term we used the fact that n2jC2j_iC2i-iP{i) = c\j_^C2i-in2jP(i) — 
because P(i) contains the term (1 — n2j)- Similarly we have n2j-iP{i) = 0. Finally, we note 

that P(^j){cljC2j-i + clj_^C2j)P{j)\0)j = and P{j){cljC2j-i + 4i-iC2i)Aj> = ^o we 
can simply replace this term by nj. 

Adding the potential and kinetic terms we obtain 

L+ L+ 

^+ = E E 2A.->4^iA^> + E 2A.> = 2//ehain. (A.4) 

j=l i=j±l j=l 

Appendix B: Numerical data on ground state degeneracy 
1. Size 3 X L, L odd 

We show the ground state degeneracy for the 3 leg ladder of odd length L for the non- 
trivial sectors with particle number / = (L±l)/2 and momenta p^; = 27r— = 27rfc/3, with 
A; = 0, 1, 2. We also indicate the total number of zero energy states, N^g, and the Witten 
index, W. We present data for L = 3, . . . , 21. 



3x3 
= 2,W = -2 



3x5 



N^, = 1,W = 1 



3x7 
N„, = 1,W = 1 
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3 X 11 


3 x 13 


3 x 15 


3 X 17 
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W = 1 


Ngs = 5, 


W = 1 


iVgs = 24,W = -2 


iVg, = 15, = 1 


{f,k) 


12 


{f,k) 


12 


{f,k) 


12 


{f,k) 


12 


5 
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7 


3 5 5 
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8 
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3 X 19 

iV,, = 23, = 1 



3 X 21 
iV^s = 132, W =-2 



{f,k) 


1 2 







1 2 


9 


11 


10 


21 


22 22 


10 


12 


11 


21 


23 23 



2. Size 3 X L, L even 



For the 3 leg ladder the total number of ground states is understood analytically as 



explained in section [III B[ The degeneracy per momentum sector can also be understood 
from the mapping between cohomology elements to tilings up to the mismatch given in 



(III. 14). For this reason we do not give all the numerical results here, instead we explain 
the momentum assignment derived from the tilings for one example. We consider the 3 leg 
ladder of length L = 12. The numerical computation of the ground state degeneracy per 
momentum sector reveals the following pattern. 



3 X 12 
N^, = 130, W = 130 





1 2 





4 3 3 
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3 4 4 
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3 4 3 


3 


6 4 4 


4 


3 3 3 
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3 4 4 


6 


4 3 3 
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3 4 4 
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3 3 3 


9 


6 4 4 


10 


3 3 4 


11 


3 4 4 



We now reconstruct this pattern via the tiling correspondence. One example of a tiling is 
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given figure III.6[ Let us refer to tliis tiling as RRGGRG, wliicli stand for the sequence of red 



(R) and green (G) tiles. The unit cell of this tiling contains 36 sites and thus corresponds to 
36 ground states. We write the eigenvalues of translations along u = (0, M) and v = {N, 0) 
as ty = e^^'^yl^^^ and = e^'^*'^-/^, respectively. It follows that the tiling RRGGRG has 
eigenvalues t^^ = 1 and ty = 1. Therefore the 36 ground states are distributed evenly over 
the 36 momentum sectors. This is also true for the ground states corresponding to the 
tilings RRRGGG and GGRRGR. We thus have so far accounted for 3 ground states in each 
momentum sector. The tiling with all red tiles and the tiling with all green tiles each have 
a unit cell of 6 sites. They each correspond to 6 ground states and taking into account the 
fermionic character of the particles we find the following for the translation eigenvalues. The 
red tiling has eigenvalues t^ty = — 1 and ty = l and the green tiling has = — 1 and ty = 1. 
The momentum assigment is thus {kx,ky) = (3, 0), (9, 0), (1, 1), (7, 1), (5, 2), (11, 2) for the 
red tiling and {k^, ky) = (3, k), (9, k) with A; = 0, 1, 2 for the green tiling. Finally, there is 
a tiling which has red and green tiles alternating leading to a unit cell, RG, of 12 sites. 
The translation eigenvalues are f^ty = 1 and ty = 1, leading to the momentum assignment: 
{kx, ky) = (3n — k mod 12, k) with A; = 0, 1, 2 and n = 0, 3, 6, 9. We readily check that this 
argument correctly reproduces the degeneracies in all momentum sectors with the exception 
of the sectors = 2tt — py = 27rA;/3 with k = 1,2. In these two sectors the degeneracy is 
3, whereas the tiling argument gives 4. It follows that of the 12 tilings of type RG only 10 
correspond to ground states. 

3. Size 4 X L 

We show the ground state degeneracy for the 4 leg ladder of length L for the non-trivial 
sectors with particle number / and momenta {Px,Py) = (27rA;^/L, 2'7r/Cj^/4). We also indicate 
the total number of zero energy states, Ngs, and the Witten index, W. We present data for 
L = 4 12. 



45 



4x4 
N^^ = 2S,W = -23 
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AT =ll,W = 11 
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TV^s = 29, W = 25 
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A/gs = 193, W = 193 
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TVgs = 151, W = -29 
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4 X 10 
293, W -- 
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4 X 11 
Afgs = 859, W = 859 



4 X 12 

TV^s = 1439, W = -1295 
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4. Size 5 X L 

We show the ground state degeneracy for the 5 leg ladder of length L for the non-trivial 

sectors with particle number / and momenta {px,Py) = {2'iTkx/ L,27iky/5). We also indicate 
the total number of zero energy states, Ngs, and the Witten index, W. We present data for 
5,6,7. 
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5x5 
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5x6 
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TVgs = 215, W = 211 
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Size 6 X L 



We show the ground state degeneracy for the 6 leg ladder of length L for the non-trivial 
sectors with particle number / and momenta {px,Py) = {27Tkx/ L,2nky/6). We also indicate 
the total number of zero energy states, Ngs, and the Witten index, W. We present data for 
L = 6,7,8,9. 
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6x8 

A^gs = 619, W = -415 
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6x9 

iVgs = 1926, W = 1462 
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6. Size 7 X L 

We show the ground state degeneracy for the system of size 7 x 7 for the non-trivial 
sectors with particle number / and momenta {px,Py) — i^-Kky/l). We also indicate 

the total number of zero energy states, A^gg, and the Witten index, W. 
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